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We present a detailed description of the nature of the wavefunction and spin distribution of the 
zero energy modes of zigzag graphene ribbons (ZGRs) in the presence of the intrinsic spin-orbit (I- 
SO) interaction. These states characterize the quantum spin Hall (QSH) phase in graphene ribbons. 
We provide analytic expressions for wavefunctions and show how these evolve as the strength of the 
interaction and the ribbon width are changed. For odd-width ribbons, we show that its insulating 
nature precludes the existence of a QSH phase. For these systems the I-SO interaction is predicted 
to have a stronger effect as shown by the enhancement of the gap as the interaction strength is 
turned on. 
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I. INTRODUCTION 

Much of the interest in graphene, the mono-layer of 
carbon atoms arranged in a honeycomb lattice obtained 
for the first time few years agc^'^ lies on its potential 
applications in electronic circuitry. As methods of fabri- 
cation improve, samples with higher mobilities and bet- 
ter conducting properties are produced, as shown by 
the recent experiments where mobilities of the order of 
20 X W^crri^ /Vs were obtained^ . It is thus reasonable 
to expect that controlled design of sample sizes as well 
as of edge terminations will be possible in a near fu- 
ture. Of particular interest for device applications are 
graphene ribbons or wires where the role of confinement 
has greater influence on transport properties. Several 
theoretical studies^i^i^i^ have focused on various proper- 
ties of graphene ribbons modeled with two different edge 
terminations -and further combinations of them to rep- 
resent more complicated edges. These standard model 
terminations are known as armchair and zigzag edges. 
While the physics of armchair ribbons shows the phe- 
nomena expected from confinement of a graphene sheet 
(the band-structure of the two-dimensional material is 
fully reproduced in the limit of an infinite wide ribbon)"*, 
zigzag ribbons possess zero energy modes that are highly 
localized along the edges of the sample. These modes, 
with a topological origin, play an important role on the 
magnetic properties of the ribbon and have been the sub- 
ject of extensive analytic and numerical researcbS. In 
this paper we present a detailed analysis of the prop- 
erties of zigzag graphene ribbons (ZGRs) that includes 
the effects of the intrinsic spin-orbit (I-SO) interaction 
introduced by Kane and Mel e -" i - - for graphene. As pre- 
dicted in that work, the I-SO interaction gives origin to 
the quantum spin Hall (QSH) phase in two-dimensional 
graphene, a new state of matter characterized by an in- 
sulating bulk and zero energy chiral modes that are fully 
spin polarized. The analytic solution of a tight-binding 
model for the ribbon allows to write explicit expressions 
for the wavefunctions as functions of the I-SO interac- 
tion strength and the ribbon width. We also are able to 



show that only even-width ribbons can exhibit the QSH 
phase and show how the gap in the energy spectrum of 
odd-width ribbons changes with the I-SO coupling. 

II. MODEL OF GRAPHENE WITH I-SO 
INTERACTIONS 

The crystalline structure of 2D graphene is described 
in terms of two sublattices A and B with a lattice con- 
stant a = 2.4 A (see Fig. ([T])). The corresponding two- 
valued wavefunctions are spinors whose components de- 
scribe the occupation of each sublattice. A nearest neigh- 
bor tight-binding model with SU(2) spin symmetry pro- 
duces a Hamiltonian matrix and corresponding wavefunc- 
tions (for a given spin) of the form: 

"-(It) w 

with (j){k^,ky) = t{e^y^^/^ + 2cos^e-*'=''''/3) ^nd h = 
a\/3/2. In these expressions (j) is given by (j){kx,ky) — 
4>{kx^ ~ky) and for real values of ky, = 0*. 

For convenience we introduce a global gauge transfor- 
mation ub ub&^^^^^^ which correspond to a redefini- 
tion of as 

0(fc„fc,)=t(e^'^-''^-f 2cos^). (2) 

The transformation amounts to label all the atoms along 
each zigzag line by a unique y-coordinate. 

The eigenvalues of ([!]) are obtained in a straightfor- 
ward manner as E = ±e = zL^/cj)^ and the corresponding 
eigenvectors are given by: 

*± = [ ) e^'^^e^'^y (3) 

where a is defined hy (f> = ee*". For neutral graphene, ^+ 
(^f-) represents solutions with E > {E < 0) and refers 
to electron (hole) conduction (valence) bands. In this 
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FIG. 1: The left panel shows a 2d graphene lattice with the 
sublattice sites labeled as A and B and the zigzag ribbon 
boundary conditions (see text). The right panel shows the 
first Brillouin zone of graphene with the two Dirac points 
chosen at A" = (|j,0) and K' = (|f , f ). 




FIG. 2: Left: I-SO term with second-neighbor hopping of 
an electron in sub-lattice A, coupled to a spin on a _B site. 
The I-SO term results from adding clock- and anticlockwise 
motions and exchanging the role of the two sublattices. Right: 
Definitions of a and p. 



language, particle-hole symmetry implies that for each 
electron state with energy E — e and eigenstate a, there 
is a hole state with E = —e and eigenstate given by a-l-Tr. 
At the Dirac points (two independent degeneracy points 
in the Brillouin zone) (j) = or cf) — 0. 

To include the I-SO interaction we use its real 
space representation in second quantization that in- 
volves a spin-dependent second-neighbor hopping term: 
Hi-SO ~ ifid^k X dkj)z4s^Cj, with d^fc as first neigh- 
bor lattice vectors connecting electrons at positions i and 
k, c| the electron creation operator at site i and the 
spin along the z-direction (we have chosen the z-direction 
perpendicular to the graphene plane for convience) ^*^'^^ . 
Notice that this expression satisfies all the symmetries of 
the graphene latticei^. Reported values for the interac- 
tion strength t' obtained with various ab-initio calcula- 
tions range from 1.2K to lOmK—. 

The total Hamiltonian in reciprocal space reads: 



H = 



S7 (j> 



(4) 



where ^{kx,ky) — 2t'{shikxa — 2shi^^^ coskyh) and 
s = ± stands for spin- up/down electron. The effect of 
the I-SO interaction as it appears in Eq.Q is to intro- 
duce opposite staggered magnetic fields acting on oppo- 
site spins. Fig.([2l) shows a schematic representation of 
one possible process introduced by the I-SO coupling. 
Below we consider only spin- up electrons. Notice that to 
obtain the corresponding expression for spin-down elec- 
trons it is just enough to make the replacement t' — > —t'. 

Spin- up electrons and holes in states fc^, ky have ener- 
gies given by -E = ±e — ±\/ (/)(() + j'^ . The main effect of 
the I-SO interaction is to open a gap at the Dirac points 
with the system becoming a bulk-insulator^*^. 

As an example, electrons eigenstates {E > 0) are given 
by: = Xk e*'^»^e*'^="^ where 



cos |e'"/2 
sin §e-'"/2 



(5) 



and the angles a, (3 are defined by 7 = ecosf] and 
(j) = |(/)|e*" . The corresponding hole state is = 



Tyj, Q^^vV Q^^''^ with 77 obtained from x by the replacement 



III. 



ZIGZAG GRAPHENE RIBBONS 



To study the role of the I-SO interaction in a confined 
geometry we analyze a zigzag graphene nanoribbon, de- 
fined according to Fig. Hard- wall boundary condi- 
tions are imposed by setting ua = on the lower border 
and ub = on the upper borderi^i^ii^. For symmetry 
reasons we chose the origin of the y-axis in the center of 
the ribbon and the boundary conditions become: 



ua{v^-WI2) = {), ub{v^W/2) = Q. 



(6) 



where W is the ribbon width with a number N — W/ 6—1 
of chains inside it. The wavefunction of the ZGR is 
found as follows: Because of translation invariance along 
the X-axis, kx is a good quantum number. For a given 
fcj;, the total wavefunction must be a superposition of 
degenerate states with different ky values. In the ab- 
sence of I-SO there are only two degenerate spinors for 



a fixed value of k^. : k,, = k and k„ 



-k. Therefore the 



wavefunction is the superposition of these two spinors: 
^I* = a'if{kx, k) + b'if{kx, —k). Application of the bound- 
ary conditions given in Eq.([S]), renders b ^ —a with 



C 



where k satisfies 



sin(a/2 + ky - mr/2) \ j^^^ 
sin(— a/2 + ky — n7r/2) 



kW 



(7) 



(8) 



and C = -y=-\/|5ftfc/ sinh(5Rfc) — sin(5fc| is the nor- 
malization factor. Notice that a peculiar feature of a 
ZGR is that k can take complex values between two Dirac 
points^i^i^. Fig. ([3]) shows the conduction bands of a rib- 
bon with W = 46. 

The expression of an edge state wavefunction with 
kxtt > TT is given by 



, siiMq{y + W/2)) 
^ ~ ' -sinh(q(y~W/2)) 



(9) 
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FIG. 3: Energy bands of a zigzag ribbon with VK = 46 in 
the absence of the I-SO interaction. Each band is doubly 
degenerate due to the spin SU(2) symmetry. The zero energy 
(edge) band (red hne) is obtained when the wavenumber k in 
^ takes an imaginary value. The edge state starts at the 
point where cos(fc°a/2) = ± 2iv^2 



where ky — —iq is purely imaginary and satisfies Eq.([8]) 
with n = and a = n + iao and ag = —qW. For k^a < tt 
the edge state wavefunction is given by 

^ - rf smh{q + iTT){y + W/2) \ ^ik.x(.(.. 
^ - y sinh[(9 + tTr){y - W/2) + ittN] J 

where ky = —ir/b — iq and satisfies Eq.® with n — N 
and a = iao and ao = —qW. 

A. Zigzag ribbons witli I-SO interactions 

Because the I-SO interaction involves second-neighbor 
hoping there is a new set of boundary conditions that 
has to be satisfied: besides those given in Eq.® it is 
necessary to impose two extra conditions given by {y — 
W/2) = usiy — —W/2) ~ 0. Correspondingly, for a 
fixed value of kx , there are now four degenerate states at 
ky = {±ki;±k2) as shown in FigU] 

Note that for some energies and values of k^ , the num- 
bers ki and k2 can be complex. To guarantee the de- 
generacy of these points, the energy condition E{±ki) — 
E{±k2) is transformed into the condition 

cosfci6 + cosfc26 = 2c(l - T^) (11) 

with = l/(8t'2sin2 ^) and c = cos(fc^a/2). 

As an example let's consider the wavefunction for a 
spin- up state with k^a > tt given by the linear combina- 
tion of the four degenerate spinors: 

+a2Xfc.e^'=^^-62X-fe.e-^'=^^. (12) 
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FIG. 4: In the absence of I-SO, for a fixed value of k^, there 
are only two degenerate states ai ky — ±k. In the presence 
of the I-SO interaction there are four degenerate states at 

ky = ±fcl, ±^2. 

Nontrivial solutions of Eq. exist only if the condition 

((71 - 72)' + (01 - '/'2)(^i - h)) sm^{qW/2) 
= ((71 - 72)' + i^i - 02)(0i - 02)) 8m\kW/2) (13) 

is satisfied. In these expressions 7j — 'y{ky — kj), and 
k = (fci+fc2)/2, q = (fci— fc2)/2. The two conditions given 
by Eqs. fTT]) and uniquely determine ^1,^2 as func- 
tions of t, t' and W. The effects of the I-SO interaction 
on the band structure have been described in a previous 
worki^. There it was shown that the quasi-degeneracy of 
the zero energy mode is lifted and the dispersion around 
fcj; = TT becomes linear. This can be shown to be the exact 
solution in the limit of a semi-infinite ribbon, i.e, when 
t' /t ^ b/W. In this case, a simplified expression for the 
edge-band dispersion can be found because the conditions 
imposed by Eqs. (fTTj) and (11311 are simplified as follows: 
cosh 2qb^l+T'^ and cosh 2kb = 4:C^{1-t^)^ /{2 + t^)~1. 
The dispersion relation of the edge sates becomes: 



e = ±6t'sin(fc^a)y 1 + 16t'2sin2(fc:ra/2) (14) 

which near k^a — tt can be approximated hy e ~ ±6t'kxa 
with the velocity of right (left) movers given by hv — 
±6t'a. 

To gain better understanding of the properties of the 
wavefunction associated with this mode, we analyze the 
behavior of the coefficients appearing in Eq. as func- 
tions of the I-SO coupling strength and the ribbon width. 
The expressions for these coefficients are: 

01 =< V-k2\X-ki >< V-ki\Xk2 > sinhfciy sinhfc2W 
bi =< V-k2\Xki >< V-ki\Xk2 > sinhgVFsinhfc2T4^ 

02 =< 'n-k2\X-ki >< V-kAxkt > smhkWsinhkW 
62 =< V-ki\Xk2 >< V-ki\Xki > sinh/cM^sinh 911/(15) 
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FIG. 5: Lowest energy band of a semi-infinite zigzag graphene 
ribbon with a I-SO coupling constant — .It. The edge 
state dispersion is given exactly by Eq. (I14p . Difi'erent colors 
represent difi'erent values of ky as described in Ref.i-. 
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FIG. 6: Coefficients of the edge state wavefunction as a func- 
tion of the I-SO coupling t' for W — 3b (left panel) and as a 
function of the ribbon width W for t' — .It (right panel). 



Notice that in the limit t' 0, e'''^'' w ~2ct'^, hi -> -62 
while al « a2 and both of them go to zero as i' ^ 
as shown in FiglGl In this way, the wavefunction of a 
ZGR in the absence of I-SO interaction is recovered. On 
the other hand for a semi-infinite ribbon: ai and fei — s- 
while 02 —62 and remains finite. 

Further insight into the effect of the I-SO interaction 
is gained when the spin distribution of the edge state 
is studied. Let us take for example spin-up states with 
k^a > TT with a given state in the conduction band (elec- 
tron or right mover state) represented by 



(16) 



The corresponding valence band state (hole or left 
mover state) is given by 



XLT = 



-VB(-y) 
<^A(-y) 



(17) 



The spin-down electron states can be obtained by the 
replacement t' —t'. The corresponding wave functions 



FIG. 7: Probability distribution for spin-up right- (filled sym- 
bols) and left- mover (empty symbols) as a function of the po- 
sition across the ribbon, for I-SO coupling strengths t' = O.lt 
(circles) and t' = O.Olt (triangles). The ribbon width is 
W = 4b. 



are: 

XRl = 

for a right mover and 

XLl -- 



v>B{y) 



(18) 



(19) 



for a left mover. In Fig. ([7]) we plot the probability dis- 
tribution for a spin-up right (spin-down left) mover as a 
function of the position y across the ribbon for two values 
of the I-SO coupling strength. Notice that the predicted 
localized^" wavefunction becomes spin-polarized as as the 
I-SO interaction is turned on. 



IV. RIBBON'S WIDTH AND I-SO 
INTERACTION 

The role of the I-SO interaction becomes more pro- 
nounced for certain widths of zigzag ribbons. Ribbons 
with even widths W = 2Mb (where M is an integer) 
have an odd number N of chains inside and are metals in 
which edge states from conduction and valence bands be- 
come degenerate at k^a = it. At the same time, ribbons 
with odd widths W = (2M -I- 1)6 (with an even number 
N of chains inside) are insulators^. Thus, it is expected 
that the effects of the I-SO interaction will be more evi- 
dent for odd-width ribbons by increasing the size of the 
gap. In fact, from Eq. pT|) it is possible to show that at 
kxCL = vr: 

kyb = TT + iqb smh{qb) = t/At' (20) 
Furthermore, for N even Eq. p^ reads 
sinh^g = (cosh 2^11/ -cosh 217) /2(1 + cosh 2(7 H/) (21) 
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FIG. 8: Band gap A = 2Eo for odd width zigzag ribbons (N 
even) at k^a — n as a, function of tire I-SO coupling t' for 
iV = 2,4,8. 

which resuhs in the foUowing energy for the edge state 
at k^a = TT 

Eo = ±tcosh{kyb)/ cosh{kyW). (22) 

For smah t' the energy of this mode increases and the gap 
(A — 2Eo) between conducting and valence bands scales 
as A w t{2t' /t)^ . In FigElwe plot a typical dependence 
of the energy Eq as a function of t' for several values of 
N. 



on the wavefunction and spin-distribution properties of 
zigzag graphene nanoribbons. Hard- wall boundary con- 
ditions in the presence of the I-SO interaction introduce a 
new set of equations to be satisfied that affect both sub- 
lattices. Thus, to obtain the wavefunction for the zero 
energy mode it is necessary to mix four different spinors 
in contrast with the two needed in absence of the I-SO 
interaction. An analysis of the role played by each spinor 
(at fixed values of the coupling strength) as the ribbon 
width is changed shows how the wavefunction in the semi- 
infinite ribbon limit is recovered. These results allows to 
track back the origin of the predicted spin-polarized chi- 
ral edge states. We also provide explicit expressions for 
the corresponding wavefunctions that fully describe the 
spin probability distribution and show an increasing spin 
polarization as the interaction strength increases. The 
effects are expected to be more immediate for odd-width 
zigzag ribbons where the I-SO interaction is predicted to 
enhance a nascent gap that makes these systems insula- 
tors. 
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